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1.0  Introduction; 

CLFE20  was  developed  under  the  assumption  of  a  state  of  generalized 
plane  strain.  The  most  general  displacement  field  for  such  problems  can 
be  written 

^total  =  ^  "x  "" 

V  =  V(y,z)  (lb) 

w  =  W(y,z)  (Ic) 

The  assumption  of  generalized  plane  strain  allows  the  user  to  study  the 
response  of  an  infinitely  long  body  of  arbitrary  cross-section  to 
various  loadings  with  the  restriction  that  all  quantities  except  axial 
displacement  are  independent  of  the  axial  coordinate  x.  The  coordinate 
system  used  by  CLFE20  is  shown  in  Fig.  (1.1).  The  x,  y  and  z  axes 
represent  the  global  coordinate  system  while  the  material  principal 
directions  are  represented  by  the  1,2  and.  3  axes.  The  region  of 
interest  to  be  modeled  by  finite  elements  is  then  a  typical  cross 
section  (or  a  symmetric  portion  of  it). 

The  element  used  in  CLFE2D  is  a  four  node,  linear,  general  quadri¬ 
lateral,  isoparametric  element,  A  three  node  version  of  the  element  can 
be  obtained  by  "double-numbering"  the  first  local  node  in  the  element 
connectivity  data.  The  available  element  geometries  are  shown  in 
Fig.  (1.2).  CLFE20  requires  that  the  local  node  numbering  scheme,  for  a 
given  element,  proceed  in  a  counter-clockwise  direction. 

An  in-depth  description  of  the  finite  element  used  in  this  program 
can  be  found  in  reference  [1].  The  finite  element  formulation  for  the 
generalized  plane  strain  problem  can  be  found  in  reference  [2]. 

In  addition  to  the  above  overview  of  the  program  capabilities,  this 
report  presents  the  input  requirements,  the  output  capabilities,  details 
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of  the  input  card  sequence,  a  note  on  symmetry  requirements  and  an 
example  problem. 

2.0  CLFE2D  Input  Requirements: 

In  this  section  a  description  of  each  type  of  card  image  in  the 

input  file  is  given.  The  required  format  of  each  card  image  is  given  in 
section  3.0. 

Type  1  Title  (Cards  1  i  2) 

The  title  consists  of  two  80  character  lines  which  will  be 
printed  on  the  output  file. 

Type  2  Master  control  Information  (Card  3) 

On  this  card  the  user  specifies  the  number  of  elements  and 
nodes  in  the  mesh,  the  number  of  different  angles  of  material 
principal  orientation,  the  number  of  specified  nodal  degrees  of 
freedom  and  forces.  The  user  also  specifies  the  desired  plot 
option  and  a  data  check  option  on  this  card  (see  section 
3.0).  The  data  check  option  permits  the  user  to  check  the 
input  without  proceeding  to  the  solution. 

3  Nodal  Coordinate  Scale  Factors  (Card  4) 

To  ease  input  the  user  can  specify  scale  factors  which 
independently  scale  the  magnitude  of  the  Y  and  Z  nodal 
coordinate  values. 

Type  4  Pri^f.-r  Control  Keys  (Card  5) 

On  this  card  the  user  sets  the  various  output  keys  (see  section 
3.0)  which  allow  for  the  generation  of  a  tailored  output  file. 

5  Plotter  Control  Information  (Card  6) 


On  this  card  the  user  specifies  the  plotter  scale  factors  (see 
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Output  Capabilities).  Card  6  can  be  omitted  if  output  plots  are 
not  desired. 

Type  6  Material  Property  Information  (Card  7  and  8) 

On  these  two  cards  the  user  inputs  one  set  of  three  dimensional 
orthotropic  material  properties.  These  properties  must  include 
thermal  or  hygroscopic  expansion  coefficients  if  the  response 
to  hygrothermal  loading  is  desired. 

Type  7  Material  Property  Orientation  Information  (Card  9) 

The  user  inputs  an  angle  number  and  magnitude  (degrees)  for 
each  different  material  property  orientation  in  the  laminate. 

Type  8  Element  Data  Information  (Card  10) 

The  element  connectivity  data,  material  property  orientation 
angle  number,  and  the  output  points  (see  section  3.0)  for  each 
element  are  specified  on  this  card.  Card  10  is  repeated  NEM 
times,  where  NEM  is  the  number  of  elements  in  the  mesh. 

Type  9  Nodal  Data  Information  (Card  11) 

The  nodal  coordinates  are  specified  on  this  card.  The  nodal 
coordinate  values  must  be  greater  than  or  equal  to  zero.  Card 
11  is  repeated  NODS  times,  where  NODS  is  the  number  of  nodes  in 
the  mesh. 

Type  10  Specified  Nodal  Displacement  Information  (Card  12) 

On  this  card  the  user  identifies  the  node  number  and  displace¬ 
ment  of  each  constrained  degree  of  freedom.  Card  12  is 
repeated  for  each  degree  of  freedom  that  has  a  specified  dis¬ 
placement.  If  there  are  no  constrained  degrees  of  freedom, 

Card  12  is  omitted  from  the  input  file. 


-  ✓'  * 


Type  11  Specified  Nodal  Force  Information  (Card  13) 

On  this  card  the  user  identifies  the  node  number  and  the  compo¬ 
nent  magnitude  of  each  of  the  specified  nodal  forces.  Card  13 
is  repeated  for  each  component  of  the  specified  nodal  forces. 

If  there  are  no  specified  nodal  forces.  Card  13  is  omitted  from 
the  input  file. 

Type  12  Normal  Strain  and  Hygrothermal  Load  Information  (Card  14) 

Here  the  user  specifies  the  applied  normal  strain  e°  along  with 
the  temperature  or  moisture  change.  It  is  possible  to  specify 
these  values  as  zero. 

3.0  CLFE20  Output  Capabilities; 

CLFE20  solves  for  the  nodal  displacements,  strains,  stresses  and 
strain  energy  densities  for  each  element  in  the  mesh.  These  values  are 
calculated  in  the  global  (x,  y,  z)  coordinate  system.  For  a  given 
element,  the  values  of  strain,  stress  and  strain  energy  density  are 
calculated  at  a  combination  of  any  nine  output  points.  These  output 
points  are  specified  for  each  element  on  Card  10  of  the  input  file. 

Figure  (1.3)  illustrates  the  location  of  these  output  points  on  a  given 
element.  Output  point  5  is  the  only  valid  output  point  for  a 
triangular  element.  Specifying  other  values  for  triangular  elements  can 
cause  errors. 

It  should  be  noted  that  the  strains,  stresses,  and  strain  energies 
that  are  printed  are  the  response  to  the  combined  mechanical  and 
Hygrothermal  loadings.  The  u  displacement  that  is  printed  is  the  U(y,z) 
of  egn.  (la).  To  obtain  the  total  u  displacement  at  a  desired  x 
location,  one  must  take  into  account  the  applied  normal  strain  e°. 
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In  CLFE2D  tne  user  has  the  ability  to  tailor  the  output  file.  By 
setting  various  printer  control  keys  on  Card  5  of  the  input  file,  he  can 
opt  to  have  all  or  any  of  the  following  printed: 


Key  (1) 

Element  Data 

Key  (2) 

Nodal  Data 

Key  (3) 

Specified  Nodal  Displacements 

Key  (4) 

Specified  Nodal  Forces 

Key  (5) 

Displacements 

Key  (6) 

Strains,  Stresses  and  Strain  Energy  Densities 

The  user  can  also  opt  to  plot  the  undeformed  and  deformed  meshes  by 
selecting  the  appropriate  plot  option  on  Card  3.  The  user  controls  the 
scale  of  these  plots  by  specifying  various  plot  scale  factors  on  Card 
6.  The  various  plot  options  and  scale  factors  are  described  in  detail 
in  the  following  outline. 

1,  Plot  Option  NPLOT  (Card  3) 

NPLOT  is  one  of  the  variables  which  is  user  specified  on  Card  3. 

The  user  determines  which  plots  will  be  printed  by  setting  NPLOT  to 
any  of  the  following  values. 

0  =  No  plots 

1  a  Undeformed  plot  only 

2  *  Deformed  plot  only 

3  «  Both  undeformed  and  deformed  plots 

2.  Plotter  Control  Information  (Card  6) 

On  Card  6  of  the  input  file,  the  user  specifies  the  following  plot 


scale  factors: 


PYSCL  =  Plot  Y-di recti  on  scale  factor 

This  variable  determines  the  maximum  length  of  the  plot  in  the 
Y-di recti  on.  PYSCL  must  be  chosen  such  that  the  maximum  Y- 
coordinate  In  the  mesh  will  be  scaled  to  a  value  less  than  or 
equal  to  18  inches. 

PZSCL  =  Plot  Z-direction  scale  factor 

This  variable  determines  the  maximum  length  of  the  plot  in  the 
Z-direction.  PZSCL  must  be  chosen  sue*'  that  the  maximum  Z- 
coordinate  in  the  mesh  will  be  scaled  to  a  value  less  than  or 
equal  to  9  inches. 

VMAX  =  Maximum  allowable  v-displacement 

This  variable  determines  the  size  of  the  maximum  v-displacement 
in  inches.  All  of  the  smaller  v-displacements  are  scaled 
proportionally  to  this  value. 

UMAX  =  Maximum  allowable  w-displacement 

This  variable  determines  the  size  of  the  maximum  w-displacement 
in  inches.  All  of  the  smaller  w-displacements  are  scaled 
proportionally  to  this  value. 


4,0  CLFE2D  Input  Card  Sequence 


Cards  IS  2  {20A4) 
Column 
1-80 


Title  Cards 
Contents 
Title 


Card  3  (1615) 
Column 


1-5 

NEM 

6-10 

NODS 

11-15 

HANG 

16-20 

NSDF 

21-25 

NS8F 

26-30 

NPLOT 
-  0 
=  1 
=  2 
»  3 

41-45 

NCHECK 

Master  Control  Card 
Contents 

Number  of  Elements 

Number  of  Node  Points 

Number  of  Different  Angles 

N'jmber  of  Specified  Nodal  Displacements 

Number  of  Specified  Nodal  Forces 

Plot  Option 

none 

undrformed  plot  only 
deformed  plot  cniy 
both  undeformed  and  deformed 
Data  Check  Option 


By  setting  NPLOT  .EQ.l  and  NCHECK  .NE.O  a  data  check  and  an 
undeformed  plot  will  be  generated.  To  obtain  the  complete 
solution,  another  run  must  be  made  with  NCHECK  .EQ.O. 


Card  4  (8D10.5) 

Column 

1-10 

SCAY 

11-20 

SCAZ 

Card  5  (1615) 

Column 

1-5 

KEY(l) 

6-10 

KEY(2) 

11-15 

KEY(3) 

16-20 

KEY(4) 

21-25 

KEY(5) 

26-30 

KEY(6) 

Scale  Factors 
Contents 

Y-Nodal  Coordinate  Scale  '^actor 
Z-Nodal  Coordinate  Scale  Factor 

Printer  Control  Card 
Contents 

Key  for  Printing  Element  Data 

Key  for  Printing  Nodal  Data 

Key  for  Printing  Specified  Nodal 

Displacements 

Key  for  Printing  Specified  Nodal  Forces 

Key  for  Printing  Displacements 

Key  for  Printing  Strains,  Stresses  and 


Strain  Energies  Per  Unit  Volume 
(if  KEY(i)  .NE.  0  -  Print) 


(NPLOT  ,£Q.  0  skip  Card  6  and  go  to  Card  7) 


Card  6  (8E10.5) 

Co 1 umn 

1-10  PYSCL 

11-20  PZSCL 

21-30  VMAX 

31-40  WMAX 


Plotter  Control  Card 
Contents 

Plot  V-Scale  Factor 
Plot  Z-Scale  Factor 
Maximum  Allowable  V-Displacement 
Maximum  Allowable  W-Oisplacement 
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Card  7  (8D10.5) 

Material  Properties  Card  A 

Column 

Contents 

1-10 

PROP(l) 

E-11 

11-20 

PR0P(2) 

E-22 

21-30 

PROP (3) 

E-33 

31-40 

PR0P(4) 

G-23 

41-50 

PROP (5) 

G-13 

51-60 

PR0P(6) 

G-12 

Card  8  (8010.5) 

Material  Properties  Card  8 

Column 

Contents 

1-10 

Prop (7) 

NU-23 

11-20 

Prop(8) 

NU-13 

21-30 

Prop(9) 

NU-12 

31-40 

Prop(lO) 

Alpha-11 

41-50 

Prop(ll) 

Alpha-22 

51-60 

Prop (12) 

Alpha-33 

Card  9  (8010.5) 

Angle  Oata  Card 

Column 

Contents 

1-10 

Ang(l) 

Angle  No.  1  (in  Degrees) 

11-20 

• 

Ang(2) 

• 

Angle  No.  2  (in  Degrees) 

m 

• 

• 

• 

• 

Ang(NANG) 

• 

• 

Angle  No.  NANG  (in  Degrees) 

Card  10  (5X.6I5) 

Element  Data  Card(s) 

Column 

Contents 

1-5 

BlanK 

(Element  Numbers  may  be  ins( 

Reference) 

6-10 

NOD(N.l) 

Node  #1  of  Element  N 

11-15 

N0D(N,2) 

Node  #2  of  Element  N 

16-20 

N00(N,3) 

Node  #3  of  Element  N 

21-25 

N0D(N,4) 

Node  #4  of  Element  N 

26-30 

lANG(N) 

Angle  Number  of  Element  N 

31-35 

ISTRS(N,1) 

Stress  and  Strain  Output  Co( 

=  0  -  None 

=  1  -  Side 

1 

=  2  -  Side 

2 

=  3  -  Side 

3 

=  4  -  Side 

4 

5 

6 

7 

8 
9 


Center 
"Node  1" 
“Node  2" 
“Node  3“ 
"Node  4” 


36-40 

41-45 


ISTRS(N,2 

ISTRS(N,3 
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71-75  ISTRS(N,9) 

’^Repeat  Card  10  NEM  Times 


Card 

11  {5X.2010.5) 

Nodal  Data  Card(s) 

Column 

Contents 

Blank 

(Node  Numbers  may  be  inserted  for  Reference) 

6-15 

Y(N) 

Y-Coordinate  of  Node  N  (This  value  must 

be  >  zero) 

16-25 

Z(N) 

Z-Coordinate  of  Node  N  (This  value  must 

be  >  zero) 

♦Repeat  Card  11  NCOS  Times  # 

#If  NSOF 

.EQ.  0  -  Skip  Card  12  and  Go  on  to  Card  13 

Card 

12  (215,010.5) 

Specified  Nodal  Oegree  of  Freedom  Card(s) 

Column 

Contents 

NO 

Node  Number  of  Specified  0.0. F. 

6-10 

lOR 

Oirection  of  Specified  0.0. F. 

(1  =  u,  2  =  V,  and  3  =  w) 

11-20 

VBDF(N) 

Specified  Displacement  Value 

♦Repeat  Card  12  NSOF  Times 

♦  If  NSBF 

.EQ.  0  -  Skip  Card  13  and  Go  on  to  Card  14 

Card 

13  (215,010.5) 

Specified  Nodal  Force  Card(s) 

Column 

Contents 

1-5 

NO 

Node  Number  of  Specified  Force 

6-10 

I  OR 

Oirection  of  Specified  Force 
(1  =  u,  2  =  V,  3  =  w) 

11-20 

VBSF(N) 

Specified  Boundary  Force  Value 

♦Repeat  Card  13  NSBF  Times 

Card 

14  (8010.5) 

Normal  Strain  and  Hygrothermal  Load  Card 

Column 

Contents 

l-lO 

EPSX 

Applied  Normal,  (X-Oi rection).  Strain 

11-20 

TEMP 

Temperature  Change  or  Moisture  Concentration 

-End- 
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5.0  A  Note  on  Symmetry  Requirements 

To  reduce  the  cost  of  a  finite  element  stress  analysis  symmetry 
arguments  are  often  used  when  building  the  finite  element  mesh.  To 
determine  the  response  of  a  thin  composite  laminate  tensile  specimen  the 
analyst  need  only  model  one  quarter  of  the  specimen  cross-section.  To 
obtain  valid  results,  however,  the  appropriate  symmetry  boundary  condi¬ 
tions  must  be  specified.  These  boundary  conditions  are  shown  in 
Fig.  (1.4).  Along  the  z-axis  (y=0)  the  U  and  V  nodal  displacements  must 
be  specified  to  be  zero.  Similarly,  along  the  y-axis  (z=0)  the  W  dis¬ 
placements  must  be  set  equal  to  zero. 


ORIGINAL  PAGE  t9 
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The  upper  right  quadrant  of  the  specimen  cross-section  will 
be  modeled  using  CLFE20 


f,2)  =  oj 


The  appropriate  boundary  conditions  for  the  quarter 
symmetry  model  are  snown  above 


Fig.  1.4  Appro pri at e_  Boundary  Conditions  for  a  Quarter 
Synwetrv  *^0061 


6.0  Example  Problem 


Consider  an  eight  ply  T300/5208  Graphite  Epoxy  laminate  with  a 
[+452/-452]s  stacking  sequence.  The  laminate  is  0.50  Inches  wide  and 
each  ply  is  0.005  inches  thick.  It  is  desired  to  study  the  response  of 
the  laminate  to  an  applied  tensile  normal  strain  =  0.001.  Of  par¬ 
ticular  interest  are  the  stress,  strain  and  strain  energy  density  dis¬ 
tributions  along  the  +45‘’/-45®  interface. 

Solution:  A  36  element,  42  node  quarter  symmetry  mesh  was  developed  to 
study  the  given  problem  using  CLFE2D.  The  input  file  and  output  for 
this  problem  are  shown  in  Appendix  A.  Note  that  the  stress  strain 
output  code  varies  for  each  element  in  the  mesh.  By  doing  this,  the 
overall  response  of  the  laminate  and  the  response  at  the  +45°,  -45° 

Interface  can  be  obtained  with  a  minimum  of  Of  put. 

A  simple  finite  element  mesh  was  used  in  the  solution  of  the 
example  problem  for  the  purpose  of  exnosition.  However,  CLFE2D  is 
capable  of  solving  large  scale  refvrch  problems. 
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ORfGJfVAL  PAGE  f3 
EXAMPLE  INPUT  DATA  QUALITY 


EXAMPLE  PROBLEM.  T300/5208.  8-LAYER, ( 45 . 0 , 45 . 0. -45 . 0 , -45 . 0 ) S 
EPSX=  0.001.  KYGROTHERKAL  LOAD  =  0.000 
36  42  4  16  0  3  0  0  0  0 

0.01  2.50D-03 

1  11  1  1  1 

32.0  100.0  0.5  0.5 

.1920D-08  1.560D+06  1.56CD-06  .4870D+06  .82000*06  .82000*06 

0.490  0.238  0.238  0.00  .1380D-04  . 1380D-04 

45.0  45.0  -45.0  -45.0 

1  4  9  10  5  1  4 

2  3  89424617 

3  2  78334938 

4  1  6  7  2  4  4 

5  9  14  15  10  1 

6  8  13  14  9  2  6  1  7 

•  7  7  12  13  8  3  9  3  8 

8  6  11  12  7  4 

9  15  23  24  15  1  5 

10  14  23  15  14  1  5 

11  14  22  23  14  1  5 

12  14  21  22  14  2  5 

13  13  21  14  13  2  5 

14  13  20  21  13  2  5 

15  13  19  20  13  3  5 

16  12  19  13  12  3  5 

17  12  18  19  12  3  5 

18  12  17  18  12  4  5 

19  11  17  12  11  4  5 

20  11  16  17  11  4  5 

21  23  32  33  24  1 

22  22  31  32  23  1 

23  21  30  31  22  2 

24  20  29  30  21  2  6  1  7 

25  19  28  29  20  3  9  3  8 

26  18  27  23  19  3 

27  17  26  27  38  4 

28  16  25  26  17  4 

29  32  41  42  33  1  5  2 

30  31  40  41  32  15 

31  30  39  40  31  2  5 

32  29  38  39  30  2  5  6  1  7 

33  28  37  33  29  3  5  9  3  3 

34  27  36  37  28  3  5 

35  26  35  36  27  4  5 

36  25  34  35  26  4  5 

1  0.0  0.0 

2  0.0  2.0 

3  0.0  4,0 

4  0.0  6.0 

5  0.0  S.O 

6  8.0  0.0 


•) 


7 

8.0 

2.0 

8 

8.0 

4.0 

9 

8.0 

6.0 

10 

8.0 

8.0 

11 

16.0 

0.0 

12 

16.0 

2.0 

13 

16.0 

4.0 

14 

16.0 

6.0 

15 

16.0 

8.0 

16 

20.0 

0.0 

17 

20.0 

1.0 

18 

20.0 

2.0 

19 

20.0 

3.0 

20 

20.0 

4.0 

21 

20.0 

5.0 

22 

2C.0 

6.0 

23 

2C  .  0 

7.0 

24 

20.0 

8.0 

25 

23.0 

0.0 

26 

23.0 

1.0 

27 

23.0 

2.0 

28 

23.0 

3.0 

29 

23.0 

4.0 

30 

23.0 

5.0 

31 

23.0 

6.0 

32 

23.0 

7.0 

33 

23.0 

8.0 

34 

25.0 

0.0 

35 

25.0 

1.0 

36 

25.0 

2.0 

37 

25.0 

3.0 

38 

25.0 

4.0 

39 

25.0 

5.0 

40 

25.0 

6.0 

41 

25.0 

7.0 

42 

25.0 

8.0 

1 

1 

0.0 

2 

1 

0.0 

3 

1 

0.0 

4 

1 

0.0 

5 

1 

0.0 

1 

2 

0.0 

2 

2 

0.0 

3 

2 

0.0 

4 

2 

0.0 

5 

■? 

4a 

0.0 

1 

3 

0.0 

6 

3 

0.0 

11 

3 

0.0 

16 

3 

0.0 

c 

3 

0.0 

34 

3 

0.0 

0.001 
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